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Abstract 

Contrary to the quark mixing matrix, the lepton mixing matrix could be symmet- 
ric. We study the phenomenological consequences of this possibility. In particular, 
we find that symmetry would imply that |C/e3| is larger than 0.16, i.e., above its 
current 2a limit. The other mixing angles are also constrained and CP violating 
effects in neutrino oscillations are suppressed, even though \Ue2,\ is sizable. Maximal 
atmospheric mixing is only allowed if the other observables are outside their current 
3cr ranges, and sin^ ^23 lies typically below 0.5. The Majorana phases are not af- 
fected, but the implied values of the solar neutrino mixing angle have some effect on 
the predictions for neutrinoless double beta decay. We further discuss some formal 
properties of a symmetric mixing matrix. 
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Low energy neutrino physics [T] is described by the neutrino mass matrix 



(1) 



where U is the leptonic mixing, or Pontecorvo-Maki-Nakagawa-Sakata (PMNS) [2], matrix 
in the basis in which the charged lepton mass matrix is real and diagonal. The three 
neutrino masses are contained in m^''^^ = diag(mi, m2, ma). A useful parameterization for 
the unitary PMNS matrix is 
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— K^vjawtj, Sij = sin 9ij and introduced the Dirac 
CP-violating phase 6. There are two independent Majorana CP-violating phases a and f3 
[3] contained in the diagonal phase matrix P = diag(l, e*°, e**-'^'^^''). Various experiments 
and their analyzes revealed the following allowed 2, 3 and Aa ranges of the mixing angles 
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Note that the present best-fit value for sin^ ^13 is zero, and that there is no information on 
any of the phases. 



Obviously, the precise form of the PMNS matrix will shed some light on the underlying 
theory of lepton flavor. One interesting possible property of the PMNS matrix is that it 
might be symmetric^, U = U^. In this letter we study in detail the consequences of this 
possibility. Recently, a symmetric PMNS matrix has been shown to follow from certain 
classes of models in which the very same unitary matrix is associated with the diagonal- 
ization of all fermion mass matrices [6]. To conduct a more detailed phenomenological 
analysis of a symmetric PMNS matrix than the one performed in Ref. [6] is one of the 
motivations of this letter. However, to put the discussion on a broader basis, let us first 
comment on the formal properties of a symmetric PMNS matrix, which are similar to the 
properties of a symmetric CKM matrixo [H [HI [9] : 

-'^Symmetry of the PMNS matrix around its Ue3-Ufj,2-Uri-8Jds has been studied in [51. 
^This possibihty has been ruled out, see below. 
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(i) first recall that in general U = Uj Uu holds, where Un is associated with the di- 
agonalization of the charged lepton mass matrix via m£m\ = Utimf'^yul with 
^diag _ (^jji^^:iji^,mr). Hence, if 

Ui = SU:, (4) 

where S is symmetric and unitary, then U is also symmetric [8]. It holds that 
5" = UiUj. If S* = 1 and mc is hermitian we recover the model from Ref. In 
this scenario we have rrii, = m^'^^ [/J and mi = U* m^^^^ [/J. Another special 
case occurs if 5* = 1 and rrif is symmetric. Hereby we obtain m* = U* mf,^^^ Ul and 
mi = U* m^'^^ ul, i.e., m* and mi are diagonalized by the same matrix; 

(ii) another formal aspect is the following [7j: we can write any unitary matrix, in parti- 
cular the PMNS matrix, as f/ = X W^^^^ X'^ , where X is unitary and f/'^'^s contains 
the eigenvalues of U: f/'^'^^s = diag(e*'^^, e'^'^^, e**^^). At first we assume non-degenerate 
eigenvalues. If X is real, then U is obviously symmetric. To turn the argument 
around, note that from U = X f/'^'^^s one can obtain U X = X [/'^"^s and - in case 
of a symmetric U - that U X* = X* W^^^^. Thus, the columns of X and X''' are eigen- 
vectors of U with identical eigenvalues. Consequently, they only differ by a phase and 
we can write X = O Q, where O is real and orthogonal and Q = diag(e*'^^, e*'^^, e*"^^). 
However, from U = X f/'^'^^s it follows that multiplying X with will lead to the 
same U and hence the phases in Q are unphysical. Thus, we have shown that a sym- 
metric PMNS matrix with non-degenerate eigenvalues implies that its eigenvectors 
are real [7]. Suppose now that two of the eigenvalues of U are degenerate: in this 
case, without loss of generality, f/'^'^^s = diag(e*'^^, e'"'^^, e*"^^), which can be written as 
f/diag ^ gi^i ^-^ ^ diagfO, 0, e*('^3-</'i) _ i)). Simply evaluating U = X U^'^'s X^ shows 

that \U\ is symmetric^. We will show next that, if \U\ is symmetric, rephasing of the 
lepton fields allows to make U symmetric; 

(iii) it is a special feature of three fermion generations that a mixing matrix having sym- 
metric moduli (lUesl = \Uri\, \Ue2\ = \Ufj_i\ and \U^s\ = \Ur2\) can be rephased in a 
way such that aigiUes) = arg([/^i), arg([/e2) = arg(f/^i) and arg([/^3) = arg(f/^2) 
[7]. To see this, consider a rephasing of the neutrino and charged lepton fields via 
z/j z/j e*°"* and £j — ^ ij e^^^ with i,j = e,fi,T or 1, 2, 3. As a consequence, the 
PMNS matrix element Uij is changed to ?7jj e**^'^*~^J^ and in addition the Majorana 
phases are modified: diag(l, e*", e**^^"'"'')) diag(e*'^^, e*^"^'^^^ g«(/3+<5+o-3)^_ Suppose the 
arguments of Uij before rephasing are In order to have arg([/jj) = aig{Uji) after 
rephasing, the parameters with which we rephase the lepton fields have to submit 
to the condition (pij — (pji = (Ji — aj + Pi — Pj mod(27r). There is a solution for this 
condition if 

Im{f/e2 U:, U:, f/.i U;,} = 0. (5) 

It is trivial to show that this equation is automatically fulfilled due to unitarity of 
the PMNS matrix in the case of symmetric moduli: consider the unitarity relation 

•^The case of all three eigenvalues being identical is trivial. 
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for the second and third row and column of U: Ue2 f^es + U*^ + Ur2 U*^ = 
and U^i U*i + U^2 Ur2 + ^^ts ^rs = 0- Multiplying the first expression with U*2 and 
the second with U*^^ and subtracting the two resulting equations while assuming 
\U^3\ = \Ur2\ yields Ue2 Ue3 U*2 = ^^/^i ^ri ^^3' This is again the condition in Eq. ([5]). 
One can show that for more than three fermion generations symmetric moduli are not 
automatically equivalent to a complete symmetric U [7]. This would be of importance 
if the LSND result was confirmed; 

(iv) a related question is the number of constraints the assumption of a symmetric PMNS 
matrix imposes on the observables. Even though there are in principle three symme- 
try conditions, l^/gsl = If^ril, \Ue2\ = lU^il and \U^3\ = \Ur2\, it is easy to see that as 
a consequence of unitarity 

\Ue3\^ - \Url\^ = |f/e2r - \U^l\^ = \U^3\^ - \Ur2\^ . (6) 

Therefore, only one constraint is inflicted on the neutrino mixing observables. 

Hence the message delivered by the last two pints is that in order to investigate the phe- 
nomenological consequences of a symmetric PMNS matrix it is obviously sufficient to 
consider symmetric moduli. This implies in particular that the Majorana phases are gen- 
erally not subject to any constraint. Only the Dirac phase and the three mixing angles will 
be affected. Moreover, the same result for the observables will be obtained for all three 
symmetry conditions. 

Note that a symmetric CKM matrix V is ruled out with current data. For instance, one 
finds that [lO] |Kb| = (3.82;°;^^) ■ lO'^ yet \Vtd\ = (8.28iJ:|) ■ IQ-^, where the errors are at 
the 3(7 level. In terms of the Wolfenstein parameterization [llj, one has \ Vub\ = AX^ [p^+v"^) 
and \Vtd\ = ((1 — p)^ + r^^). Since p 7^ | the two elements can not be equal. 

From now on we will focus on the phenomenological consequences of a symmetric PMNS 
matrix. Though one could use a parameterization suitable for the study of symmetric 
matrices and identify the elements of this parameterization with the usual mixing angles 
and CP phases from Eq. we will focus in this letter on the usual parameterization of 
Eq. ([2]) and directly obtain correlations between the neutrino mixing observables 612, 623, 613 
and S. First, let us obtain the ranges of the individual elements of the PMNS matrix: we 
vary the mixing angles in their allowed ranges given in Eq. ([3]) and the phase 6 between 

"^Such a parameterization will be useful if the PMNS matrix turns out to be indeed (close to) symmetric. 
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zero and 27r, resulting in 
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and |Jcp| < 0.037 (at 2a) 



and \Jcp\< 0.048 (at 3a) , (7) 



and I Jcp I < 0.056 (at 4a) 



We have also given the maximal possible value of the Jarlskog invariant Jqp to which any 
CP violating effect in neutrino oscillations is proportional [T2] : 



Jcp = Im {f/ei f^e2 ^Mi} ^ o siu 26'i2 siu 26123 sin 26^13 cos 6^13 sin 5 . 



Note that the rephasing of the PMNS matrix elements as discussed before Eq. 
Jcp invariant. 



leaves 



Consider now in Eq. ([7]) the symmetry condition {Uesl — \Uti\- Apparently, to fulfill 
this condition the current 2a ranges of the observables do not suffice. This can be easily 
understood qualitatively since |?7e3| is given by sin 6^13 and is therefore small, whereas 
I^tiI = |'5i2S23 — C12 C23 Si3 e^'^l is generally large. For a small \Ut-i\ of order |?7e3| it is 
necessary that ^13 is large and that S lies close to zero or vr in order to subtract the second 
term in Uri from the first one. Moreover, S12 S23— C12 C23 S13 is smaller when C23 is larger than 
S23, i.e., atmospheric neutrino mixing will tend to be governed by sin^^23 < 1/2. These 
statements can be made more precise: it is easy to proof that all three symmetry conditions, 
\Ue3\ = \Uti\, \Ue2\ = \U^i\ and |t/^3| = \Ur2\ fulfilled for one single condition, namely 

, sin 6*12 sin 6^23 , , 

|f^e3| = , „ . (9) 

V 1 — sin 6 cos^ 6*12 cos^ 6*23 + cos 6 cos 6^12 cos 6^23 

Interestingly, this constraint can also be derived by setting the real and imaginary parts 
of Ue3 = Uri equal. One gets then two relations, S13 = S12 523/(0^ + C5C12C23) and 
(s/3 + C12 C23 ss) S13 = 0. From the first one it follows that vanishing of S13 in not pos- 
sible, since experiments show that S12 and S23 cannot be zero. Realizing this, the second 
relation is fulfilled when sin/3 = —sin 5 cos ^12 cos 6*23. Inserting this expression in the 
first condition yields Eq. (Q. However, care has to be taken when working with real and 
imaginary parts of mixing matrix elements, since their individual phases have no physical 
meaning. 

With the 2 (3 and 4) a ranges of 612 and 6'23, and with varying the phase S between zero 
and 27r, one obtains from Eq. that {Uesl^ ~ 0-035 (0.028 and 0.023), which has to 
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be compared to the experimental upper limits of 0.025 (0.041 and 0.058). Therefore, a 
symmetric PMNS matrix predicts that \Ue^\ should be above its current 2cr limit. Conse- 
quently, the scenario is easily falsifiable, since the indicated value of |?7e3| should be verified 
in upcoming measurements, in particular by the Double Chooz experiment [13] (see also 
|14]): according to Ref. [I3], data taking can start in 2008 and the Sa limit on |?7e3| will 
be improved from its current value 0.04 to 0.01 (0.006) after 2 (6) years of data taking. 
These numbers are well below the prediction of a symmetric PMNS matrix. Fixing 
to 7r/4, gives |t/e3p ~ 0.050 (0.046 and 0.042). This means that maximal atmospheric 
neutrino mixing would be at roughly 3 standard deviations in conflict with a symmetric 
PMNS matrix. Fixing 6*12 to its best-fit point gives |t/e3p ^ 0.041 (0.036 and 0.033), which 
is hardly compatible with the current 3(T limit of |?7e3p- Moreover, maximal atmospheric 
mixing is not compatible with sin^ 612 = 0.30. We conclude that the current best-fit values 
of the oscillation parameters are not compatible with a symmetric PMNS matrix. To be 
precise, the prediction of sin^6'i2 = 0.30 and sin^6'23 = 1/2 would be |f/e3p ~ 0.059. To 
quantify the compatibility of a symmetric mixing matrix with current data we have also 
performed a simple minimization: we introduce 



where (s^j)best-fit and aij are the best-fit values and errors from Eq. (I3l). Obeying the 
symmetry condition Eq. (Q, one can find a minimum of = 10.29 for the parameters 
sin^6'i2 — 0.28, sin^ 6*23 — 0.36, |t/e3p — 0.035 and 5 ~ 0. The corresponding pulls for 
sin^6'i2, sin^^23, and |f/e3p are —0.81, —1.78 and 2.55, respectively. 

In Fig. [1] we show plots of sin^6'i2, sin^ 6*23 and \ Jcp\ against |f/e3|, obtained from Eq. ([9]). 
When we simply vary the mixing angles 6'i2, 13.23 in their allowed 3 and 4a ranges from 
Eq. ([3]) and require the PMNS matrix to be symmetric, the plots look identical. Note 
that the 3a ranges of the oscillation parameters imply sin^ ^23 < 1/2. In the plot of |?7e3| 
against Jcp, we also indicated the maximal \Jcp\ allowed by current data. This serves 
to illustrate that CP violating effects are rather small when the PMNS matrix is sym- 
metric, even though {Uesl is sizable. Indeed, the scenario is compatible (at 3a) with CP 
conservation, in which case |t/e3p ~ 0.035 (0.028 and 0.023). Note further that, for the 3a 
ranges of the oscillation parameters, sin^ 612 takes values on the lower side of its allowed 
range. This has interesting implications for the effective mass (m) = | ^ U^^ rrii \ governing 
neutrinoless double beta decay. If neutrinos enjoy an inverted ordering, the minimal value 
of the effective mass is {m)'^^ = cf^ \/ Am\ cos26'i2. Hence, the smaller 6*12, the larger the 
minimal value of (m) in the inverted ordering. This in turn simplifies distinguishing the 
normal from the inverted hierarchy with neutrinoless double beta decay, or fully probing 
the inverted ordering regime To quantify this statement, the lower limit in case of an 
inverted hierarchy is in general {tu)^^ ~ 0.2 c^g a/ Am^, where we have inserted the lowest 
possible value of cos26'i2 at 3a. With the constraint stemming from a symmetric PMNS 
matrix we see from Fig. [T]that sin^ 6*12 ^ 0.32 and consequently ("^)[^° — 0.36 cf^ a/ Am\. 
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This is almost a factor of two larger than without the constraint from symmetry. 

To summarize, we studied what consequences would arise from a symmetric PMNS matrix 
U, a scenario which in contrast to a symmetric CKM matrix is still compatible with data. 
We noted that in this case either the eigenvectors of U are real or U has two degenerate 
eigenvalues. A symmetric U arises when Ui and 11^ are connected by a symmetric and 
unitary matrix S via Ui = SU*. One simple example is when the neutrino mass matrix 
and the complex conjugate of a symmetric charged lepton mass matrix are diagonalized by 
the same matrix. Symmetry implies one constraint on the neutrino oscillation observables, 
not on the Majorana phases. The scenario is easily falsifiable in the near future, since 
it predicts that |?7e3| is larger than its current 2a limit of about 0.16. Experiments like 
Double Chooz can therefore easily rule out a symmetric PMNS matrix. In addition, there 
are interesting and testable correlations between the observables, as given in Eq. (Q and 
illustrated in Fig. [TJ If the 3a ranges of the oscillation parameters are taken, 623 cannot be 
maximal and lies below 7r/4. Solar neutrino mixing is far away from its maximal allowed 
value, which affects predictions for neutrinoless double beta decay. In general, the CP 
phase 6 is close to a CP conserving value. 
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Figure 1: Plots of the oscillation observables for a symmetric PMNS matrix. The correla- 
tions are a consequence of Eq. ([9]). We allowed the parameters to vary in their current 3 
(left plots) and 4cr (right plots) ranges, which are indicated in the plot. 



